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Abstract. Our ability to correct the observational photometry of galaxies depends upon our knowl- 
edge of the attenuation of light produced by the dust contained in the interstellar medium. We 
will present a model based on the statistical properties of the local density in isothermal turbu- 
lence which might be appropriate to calculate the radiative transport through the diffuse interstellar 
medium. The model will be applied to study the attenuation caused by a distant foreground screen 
and a non scattering slab where the sources are mixed within the turbulent medium. It will be shown 
how the turbulent structure affects the attenuation curve and how the attenuation curve varies with 
the thickness of a foreground screen. 



INTRODUCTION 

It is a well known phenomena that the collective star light from galaxies, in particular 
at UV and optical wavelengths, is obscured by dust grains distributed in the diffuse 
interstellar medium. To derive accurate measurements of important parameters such as 
the star- formation rate as a function of redshift, the famous Madau-plot HH], a correction 
can be essential but is still uncertain. 

For sake of clarity we want to differentiate between the dust extinction in the case 
of point sources such as stars and the attenuation in the case of extended sources such 
as galaxies. In this scheme the dust extinction is caused by dust absorption and dust 
scattering along the line of sight and is simply proportional to the wavelength dependent 
extinction coefficient and the column density A^. In contrast the dust attenuation is the 
wavelength dependent loss of the intrinsic light of all stars and can be explained by an 
effective optical depth. This optical depth is determined by several different conditions 
such as the relative distribution of emitting stars and the attenuating dust, the viewing 
angle to the galaxy, the variation of dust properties with radiation field and density 
throughout the ISM, and the structure of the ISM itself, which is far from homogeneous. 
In addition it might be important that light scattered out of the observed direction might 
be partly compensated by photons scattered into the observed direction. It is no surprise 
that, due to its complexity, the attenuation of the starlight of galaxies is still a rather 
unsolved problem. 



MODEL 



Many different attempts have been made to take the inhomogeneous nature of the ISM of 
galaxies into account either to model the spectral energy distributions 10,0] or to derive 
corrections for dust attenuation (^,'^,'6]. The model presented in the following is based 
on a realistic description of an inhomogeneous structure of the local density caused by 
turbulent motion. The basic idea of this model is the fact that the transmission of a distant 
foreground screen is, if we neglect any variations of the optical properties, primarily 
determined by the probability distribution function (PDF) of the column density. By 
using a simplified description of the turbulent density structure this PDF can be related 
in a simple way to its statistical properties and the thickness A of the dusty screen in 
front of the emitting sources. 



Model of the Isothermal Turbulent Medium 

The turbulent medium is taken to be isotropic and we assume that inside a certain 
interval the power spectrum P{k)=\p{k) p of the local density p(r) can be described by a 
simple power law P{k) <^ k" with power n. The scaling relation extends from a minimum 
scale Lmin to a maximum scale Lmax- 

As has been shown by using hydrodynamic simulations of turbulence in compressible 
fluids the PDF of the local density can be described by a log-normal density distribu- 
tion if the turbulence is approximately isothermal [1^^^. As found analytically i fioll . 
the log-normal density distribution should be the exact solution if supersonic isother- 
mal turbulence is considered. The PDF of the local density p in the ISM is therefore 
approximated by: 

pip) = lp(lnp) = ^=^e^^/< (1) 
P V27rap 

with.x: = Inp — Inpo (2) 

where Oinp is the standard deviation of the log-normal density distribution and where po 
is defined by the mean value In (p) = Inpo + 0.5a^jjp. 

In general the density contrast increases as consequence of higher compression result- 
ing from higher Mach numbers M. In particular for non-magnetised forced turbulence it 
has been found by using 3D-simulation [8, 10] that the standard deviation 

Op = (p) \/e-'^<-\ (3) 

of the local density is almost linearly correlated with the Mach number by 

ap^|3M(p) (4) 

where |3 ~ 0.5. The more general case appears to be more complicated as for magnetised 
turbulence no simple correlation between density contrast and Mach number has been 
found. As might be also expected from the additional pressure support to the plasma 
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FIGURE 1. Standard deviation of the column density, given by the ratio <^n/{n) /^p/{p)' function 
of shce thickness A/Lmax- The power spectrum of the local density is taken to be Kolmogorov with 
n = — 10/3. Also shown is the dependence on the dynamic range A;max/^min taken to be lO', 10^ 10^ 
and 10^. The long dashed line shows the behaviour in the limit of thick slices where the thickness is larger 
then the maximum turbulent scale Lmax- 



provided by magnetic fields, it seems that the density contrast becomes weaker when 
magnetised turbulence is considered LlO-.ll.l . 

The variance of the local density is related to the power spectrum of the turbulent 
density structure by: 

PDF of the Column Density 

The variance of the column densities. Equivalent to the variance of the local density 
the variance of the column density is given by the total power of the column densities: 

In general the power spectrum of the column densities along the z-axis is determined by: 

P{K) = dhP{K,k,)W{k,A) (7) 

with ^(fc^A) = I V7(fc^A) p where W{kyA) is the Fourier transform of the window function 
which defines the shape of the slice through the turbulent density structure For 
a simple box-window, as used here, we have W{k^A) = sin^(fc^A/2) / {k^A/2)-^. 

We want to call a slice through the turbulent medium thick or thin if its thickness is 
either larger or smaller than the maximum scale length Lmax- If the slice is thin with 
A <^ Lmax the power spectrum of the column density is given by P{K) oc K"+^ for 
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FIGURE 2. PDF of the column density N / (N) through a simulated density structure with a log- 
normal density distribution of the local density. The power spectrum of the local density is assumed to be 
Kolmogorov with n = —10/3 and the standard deviation of the local density Op/(p) — 1-5. The dynamic 
range is chosen to be Lmax/imin = 5.4. The slice thickness through the density structure is taken to be 



l/Lmax < ^ ^ 1/A. On the other hand, if the thickness of the slice is larger than the 
maximum 'cloud' size (A > Lmax) the power spectrum is well described by P(K) c< K" 
and therefore by the same power law as the power spectrum of the local density 1 12Lll3ll. 

In the limit of very thin slices smaller than the minimum 'cloud' size the fluctuations 
of the column densities are the ones of the local density. Due to averaging effects in 
the line of sights the fluctuations of the normalised column densities decrease with 
slice thickness (Fig.[Tl). For slices thicker than the maximum cloud size the variation as 
function of slice thickness becomes a simple power law with CJ/^/(A?)/cJp/^p) A^^/^. The 
actual value of the standard deviation of the column density depends on the statistical 
properties. In a medium with a larger dynamic range the averaging effect in the line of 
sights is larger and the fluctuations of the column densities therefore weaker. 



The functional form of the PDF. Based on MHD simulation of isothermal forced 
turbulence it has been found ifllll that the functional form of the PDF of the column 
density is approximately also log-normal. The convergence of the log-normal density 
distribution of the column density into a Gaussian density distribution in the limit of 
thick slices has been analysed in some detail by Vazquez-Semadeni and Garcia |[l4ll . 
They found by using density structures that the PDF is skewed in comparison to the 
ideal log-normal density distribution. 

A more quantitative approach has been presented by Fischera and Dopita based 
on simulated density structures of an idealised isothermal turbulent medium with a 
log-normal density distribution of the local density and where the power is given by 
a defined power law. The obtained distributions can be well approximated by a log- 
normal density distribution (Fig. El). However, deviations do exist as the fitted log-normal 
density distribution predicts too high probabilities at low but too low probabilities at high 
column densities. The deviation increases for broader distributions of the local densities. 
In addition the standard deviation of the fitted distributions is in general lower than the 
exact value fls']. Despite these uncertainties it seems to be a first good approximation 
to simplify the functional form of the PDF by a simple log-normal density distribution 
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FIGURE 3. The parameter space of the distant turbulent screen given by the mean attenuation (Ay) 
and the standard variation cTav/(Av> — '^n/{n) of the column density. The attenuation curves A^/Ay are 
well defined by the absolute to relative attenuation Ry shown as thin solid lines. The actual attenuation 
Ay is shown as the thick solid lines. 



where the variance is simply determined by the thickness of the slice and the statistical 
properties of the local density. 



We first consider a situation where the stars are seen through a turbulent dusty screen 
which is distant from the observer and the emitting stars. In this case the scattering 
makes only a small contribution to the collected star light and can be neglected. For 
simplicity we assume that the dust properties are not varying with local density or 
radiation field so that the optical depth is proportional to the column density. As the 
underlying extinction curve we use the mean extinction curve of our Galaxy which we 
adopted from Weingartner & Draine [15]. 

At each wavelength the turbulent screen leads on average to an effective extinction 
given by 



where p{y) is the normal density distribution and where e^' = x/ (x). 

In general a turbulent screen is more transparent in comparison to a homogeneous 
screen and the attenuation curve flatter than the underlying extinction curve. The effect 
of the turbulent screen on the attenuation increases with dust content or (Ay) and wider 
distributions of the column density (Fig. El). The flattening can be characterised by the 
absolute to relative attenuation Ry = Av/E(B — V) which is higher in case of flatter 
attenuation curves. In case of the foreground screen the Ry- value is furthermore a perfect 
parameter to determine the curvature of the attenuation curves (Fig.|4l) over the whole 



The Distant Foreground Screen 
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FIGURE 4. Attenuation curves A^/Ay and reddening curves E(k — V)/E(B — V) of a distant fore- 
ground screen. The curves are well defined by the absolute to relative attenuation Ry which is varied from 
3.5 to 10 in steps of ARy = 0.5. The thick line is the assumed extinction curve given by Weingartner and 
Draine 1 15|. The actual curves do not strongly depend on the absolute attenuation in V which is taken to 
be Ay = 10. The broken hne is the so called 'Calzetti-extinction-curve'. 

wavelength range lH^ .^ The dependence on the mean attenuation (Ay) on the other hand 
is only weak. In the limit of small fluctuations with <3\nAv < 1 the relative attenuation 
Ax/ {Ax) and the curvature depend only on the product of^^^ (Ay) fl^ so that in highly 
optically thick media also small fluctuations can lead to flatter attenuation curves. 



The 'Calzetti Extinction Law ' 

For the nearby universe we have some knowledge about the dust obscuration. The so 
called 'Calzetti-extinction-law' iUtI] has been obtained from a huge sample of star-burst 
galaxies and can be used to derive the intrinsic spectrum of those galaxies. The curve is 
flatter in comparison to the mean extinction curve in our galaxy and also flatter than the 
ones derived for the Large and the Small Magellanic Cloud. Li addition the feature at 
2200 A which is very prominent in case of the mean extinction curve of the Milky Way 
seems to be absent. The absence is generally explained by the destruction of the carriers 
of this feature by the higher UV-radiation field in those galaxies. The flatter curvature in 
the optical is thought to be related to the clumpy medium in front of the stars. 

The overall curvature of the Calzetti-Curve can be naturally explained by a turbulent 
foreground screen [18]. The %^-fit based on photometric measurements from star-burst 
galaxies is excellent (even though the minimum of 0.048 suggest an overestimate 
of the errors) (Fig.|5l). The fit provides a most likely Ry-value of 4.75 ±0.45 which is 
in excellent agreement with the original value Ry = 4.88 ±0.98 derived for star-bursts 
ifioll . Based on IR-measurements this value has been corrected to Ry = 4.05 ±0.80. The 
agreement with our value of a turbulent screen, however, is still good. 



In general the curvature depends also on the geometry so that attenuation curves of different geometries 
but the same Ry-value might deviate. 
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FIGURE 5. Result of the X^-fit to the 'Calzetti extinction curve' . The 'Calzetti-curve' is shown as thick 
soUd Une in the left hand figure. Also shown is the mean extinction curve of our galaxy (Weingartner and 
Draine ifisll ') used in the model. The fit has been applied to the data point which are based on Calzetti flgil 
and Leitherer et al. f^o']. The best fit is shown as dashed line. Regions of confidence are given as filled 
contours (68, 90, 95, 99%) on the right hand figure in the parameter plane defined by Ay and aAy/(Av>- 
The dashed lines show the range of Ay determined for star-burst galaxies. The dashed-dotted line is the 
mean value. 



The measured range of Ay of star-burst galaxies provides a standard deviation of the 
optical depth (column density) in the range 0.66 to 1 1 .2. If we apply the relation between 
Mach-number and standard deviation we obtain a minimum Mach number in the range 
1.3 to 22. 

It might be useful to compare those numbers with the expectation of the cold neutral 
medium (CNM) and the warm neutral medium (WNM) in our own galaxy. The velocity 
dispersion for the CNM and for the WNM has been measured (based on HI and CO 
observation, respectively) to be 6-8 km/s and 7-10 km/s. Assuming a temperature of 100 
and 6000 K for the CNM and the WNM the Mach number is roughly 12 and 1.8. If we 
use the relation El the expected range for Op/^p^ is in the range 0.9 to 6 consistent with 
the Mach numbers found by fitting the Calzetti curve by a simple turbulent screen. 



The Effect of Slice Thickness on the Attenuation Curve 

For a better understanding of a turbulent density structure on the attenuation it is 
important to know how the attenuation changes with the thickness of the dust layers 
in front of the stars. As an application one might think of a simple geometry of disk- 
like galaxies where the stars are embedded in a turbulent dusty layer and therefore 
seen through a turbulent screen. If the medium would be homogeneous the optical 
thickness of the dusty layers in front of the stars would increase with inclination angle 
as 1 / cos i. The dependence is somewhat weaker if the medium is turbulent as the 
relative attenuation \x/ i^x) decreases with the thickness of the screen. In addition 
the attenuation curves should flatten with viewing angle as the Ry-value increases as a 
function of A/Lmax- 
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FIGURE 6. Variation of the absolute to relative attenuation Ry as function of the slice thickness 
A/Ljnax for a certain choice of the turbulent density structure. The thick solid lines refer to turbulent 
media with different mean attenuation (Av), varied from 0.1 to 10. The absolute attenuation Av is 
shown as thin solid lines. The absolute Ry-values in the limit of thick slices are given as dotted lines. 
Also shown are the maximum, the mean, and the minimum attenuation Ay, derived for star-burst galaxies 
(dashed, dashed-dotted, and dashed line). 



In the limit of thick slices (or high inclination angles i) the relative attenuation 
Ax/ {Ax) and the absolute to relative attenuation Ry approach asymptotically a mini- 
mum and a maximum value, respectively. The limiting values are, for a given power 
spectrum and dynamic range, only determined by the product (xy)/^ ^o/(o) 



The Relation between Ry and Ay 

The relation between the two observable quantities Ay and Ry has consequences 
for the interpretation of the change of an intrinsic spectral energy distribution from 
galaxies caused by dust obscuration. A stellar spectrum well described over a certain 
wavelength range by a simple power law Ix °^ with power [3 i2llll will have a different 
functional form if seen through a dusty screen. If the spectrum is observed at two 
different wavelengths Xi and X2 the inferred power P would be different from P by 

A R = R _ R = '^^1 ~ ^ _}_ (Q\ 

^ ^ ^ ln?ii-lnX2 RAln?ii-lnX2' ^ 

Here, we introduced the absolute to relative attenuation R^^ = XX2/ i^h "'^A-a)- '^A.i ^nd 
XX2 the effective optical depths at wavelength Xi and X2. 

Due to turbulence the variation of A(3 is in general not directly proportional to A;^^. 
As R^^ > Rx2 the change of the spectral index P should be furthermore less strong in 
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FIGURE 7. Relative attenuation Ay / (Ay) and absolute to relative attenuation Ry as function of mean 
optical depth (x) through a turbulent slab in which the sources are homogeneously distributed. The dotted 
curves show the behaviour in case of a homogeneous medium. The solid lines correspond to different 
assumptions of the mean optical depth at one maximum scale Lmax- The labels refer to log 



comparison to a homogeneous screen. 



The Non-Scattering Turbulent Slab 

To demonstrate the effect of the geometry on the attenuation curve we consider the 
extreme case that the stars are homogeneously distributed throughout the turbulent slab. 
For simplicity the scattering is ignored so that the dust layers in front of the stars at each 
depth can be approximated by a distant foreground screen. 

The slab-geometry leads even in case of the homogeneous medium to a lower attenu- 
ation and to a flatter attenuation curve (Fig.|71). In the limit of an optically thin and opti- 
cally thick limit the effective optical depth is x = 0.5 (x) and x = In (x), respectively. The 
Ry-value increases strongly when the slab becomes optically thick. In the limit of an op- 
tically thick slab the absolute to relative attenuation varies as Ry = In (x) / ln( 1 + 1 /Rv) 
and therefore is directly proportional to the effective optical depth x. 

By introducing a turbulent density structure the medium becomes more transparent 
and the relative attenuation decreases in comparison to the homogeneous slab (Fig. |71). 
In the region (x) <~ 4 the turbulence furthermore leads to a flatter attenuation curve. 
The behaviour at higher mean optical depth depends on the dust content per maximum 
scale length. In rather optically thick environments the attenuation curve can be slightly 
steeper than the very flat attenuation curve expected for a homogeneous medium. 



CONCLUSION 

The described model relates the dust attenuation to the physical properties of the tur- 
bulent ISM and therefore might be a solution to obtain accurate corrections important 



to measure the star-formation rate against redshift. As shown, a distant turbulent screen 
can naturally explain the attenuation curve derived for star-burst galaxies and may be 
also applicable for other star-forming regions. 

The idealised model of a turbulent density structure might also be helpful to under- 
stand the radiative transport through the turbulent ISM in the disks of normal galaxies. 
The model may furthermore provide, if used in a full radiative transfer code which also 
includes the thermal emission from dust grains, a better understanding of the SED of 
star-burst galaxies or Hll-regions. 
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